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ABSTRACT 


The  application  of  queueing,  or  waiting-line,  theory  to  the  review 
and  validation  of  Foreign  Military  Sales  (FMS)  cases  by  the  Cost 
Analysis  office  of  an  Anny  Materiel  Development  and  Readiness  Command 
major  subordinate  command  Is  presented  In  the  report.  Queueing 
analysis  Is  applied  to  determine  the  behavior  characteristics  of 
the  processing  system  so  that  management  can  take  appropriate 
action  to  reduce  the  total  time  In  the  system.  The  primary  charac- 
teristics analyzed  are  the  expected  length  of  the  queue,  the  expected 
number  of  cases  In  the  system,  the  expected  waiting  time,  and  the 
expected  turnaround  time.  Queueing  theory  Is  described,  and  the 
application  to  FMS  cases  presented,  to  also  serve  as  an  aid  In  the 
solution  of  similar  types  of  waiting-line  problems. 
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A.  Purpose 


This  paper  is  presented  to  describe  the  use  of  queueing,  or 
waiting  line,  theory  in  the  analysis  of  turnaround  time  in  the 
review  and  validation  of  Foreign  Military  Sales  (FMS)  cases  within 
the  Data  Analysis  and  Validation  Branch  (DRSAR-CPE-D)  of  the  HQ, 

ARMCCM  Cost  Analysis  Division.  The  method  of  analysis  described 
is  applicable  to  a variety  of  similar  work  situations. 

B.  Queueing  Theory 

An  inportant  class  of  management  problem  includes  those  that 
can  be  characterized  as  "arrival  and  departure"  problems.  These 
types  of  problem  occur  whenever  randomly  arriving  customers  are 
required  to  wait  for  some  conmon  type  of  service.  The  term  customer 
can  denote  people,  paperwork,  machines,  or  any  other  kind  of  discrete 
arrival,  whereas  the  term  service  applies  to  any  kind  of  work  per- 
formed on  or  for  the  customer.  Exanples  are  customers  arriving  at 
a doctor's  office,  machines  waiting  to  be  repaired,  and  paperwork 
waiting  to  be  processed.  Given  that  the  service  facility  or  facilities 
are  adequate  to  meet  the  danand  for  service,  and  are  not  excessive 
relative  to  the  demand,  a condition  exists  in  which  both  customers 
and  service  facilities  occasionally  encounter  temporary  waiting  owing 
to  variations  in  service  times  and/or  arrival  rates.  IVhen  customers 
are  required  to  wait,  they  form  a waiting  line,  or  queue ; hence,  the 
analytical  techniques  used  to  solve  waiting- line  problems  are  embodied 
in  what  is  conrionly  referred  to  as  queueing  theory. 

Patterns  of  service  and/or  arrival  rates  are,  under  certain 
assunptions,  susceptible  to  approximation  by  mathematically  defined 
frequency  distributions.  On  this  basis,  mathonatical  models  have 
been  formulated  which  can  be  used  to  determine,  on  a probabilistic 
basis,  the  behavior  characteristics  of  arrival -departure  systems  that 
meet  the  required  conditions  and  assumptions.  Repeated  execution  of 
the  models,  varying  the  number  of  service  facilities  or  stations,  and, 
if  feasible,  the  mean  (average)  service  time,  provides  the  behavior 
pattern  of  the  system  being  analyzed.  These  data  can  be  used  to 
indicate  those  actions  management  can  take  to  minimize  the  total  time 
customers  spend  in  the  system,  and/or  the  total  cost  of  waiting. 

The  siirplest  kind  of  waiting  line  systan  are  those  in  which 
customers  arrive  at  a single  service  facility,  occasionally  have  to 
wait  for  service,  and  are  serviced  according  to  some  priority  rule 
such  as  "first-come,  first- served".  This  type  of  system  can  be  called 
a single- station,  single-stage  systan.  More  complex  kinds  of  system 
can  exist,  such  as  single-station,  multiple-stage;  multiple-station, 
single-stage;  and  multiple- station,  multiple- stage,  or  series  parallel. 
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The  models  formulated  for  queueing  analysis  vary  depending  on  the 
type  of  system,  although  each  has  its  basis  in  the  type  of  arrival - 
rate  and  service-rate  frequency  distributions  assumed  for  waiting 
lines.  The  type  of  model  and  tlie  input  parameters  which  are  appli- 
cable to  the  analysis  of  FMS  cases  are  described  in  paragraph  C3. 

C.  Application  to  FMS  Cases 

1.  Purpose  of  the  Analysis 

Queueing  theory  has  been  applied  to  the  review  and  valida- 
tion of  FMS  cases  by  DRSAR-CPE-D  to  provide  management  with  the 
behavior  characteristics  of  the  review  and  validation  process  so 
that  action  could  be  taken  to  reduce  the  total  expected  time  an 
individual  case  spends  in  the  system  from  receipt  to  completion, 
i.e.,  the  expected  turnaround  time.  The  relevant  characteristics 
are  the  average  time  in  the  system  and  the  number  of  reviewers 
required  to  reduce  or  minimize  the  average  time. 

2.  Type  of  System 

Cases  arrive  randomly  from  the  HQ,  AR^CCM  International 
Logistics  Directorate  (DRSAR-IL).  Following  receipt,  they  are 
logged  in  and  assigned  to  one  of  several  analysts.  A case  may  or 
may  not  have  to  wait,  depending  on  when  it  arrives  relative  to  the 
backlog  and  cases  in  process  in  the  Branch  at  the  time.  In  general, 
FMS  cases  are  reviewed  and  validated  to  ensure  conpliance  with 
pricing  regulations  and  policy,  verify  the  correct  application  of 
inflation  indices,  guarantee  full  cost  recovery,  assure  fair  prices 
to  both  parties,  ensure  proper  calculation  of  unfunded  costs,  and 
to  reduce  the  number  of  amendments  necessary  because  of  price 
fluctuations.  The  same  basic  process  is  followed  in  the  review  and 
validation  of  each  case,  although  they  differ  in  scope  and  hence  vary 
randomly  in  the  processing  time  required.  In  actual  practice,  a case 
may  be  started  and  later  delayed  pending  receipt  of  essential  infor- 
mation; hence,  an  analyst  may  have  more  than  one  case  under  review 
at  any  one  time.  Following  the  review  and  validation  process,  formu- 
lation of  the  analyst's  initial  recommendation,  and  Branch  Chief 
review  and  recommendation,  decision  for  formal  action  is  made  by  the 
Division  Chief,  Comptroller,  or  Deputy  Comptroller,  depending  on  the 
dollar  value  of  the  case  and  whether  concurrence  or  non- concurrence 
is  indicated.  The  case  is  then  released  for  return  to  DRSAR-IL. 

The  puipose  of  the  analysis  is  to  determine  the  characteristics 
of  the  analyst's  review  and  validation  process,  not  the  approval 
process.  This  restricts  the  definition  of  the  system  to  the  latter, 
and  sinplifies  the  analysis,  since  a separate  approval/decision  step 
creates  a multiple- stage  system.  As  a result,  the  output  of  the 
analysis  can  only  guide  management  actions  to  reduce  tumaround  time 
through  conpletion  by  the  analyst.  The  contribution  to  total  elapsed 


time  through  release  to  DRSAR-IL  created  by  the  approval  stage  is 
subject  to  separate  analysis.  The  review  and  validation  system 
thus  described  is  a multiple- station,  single-stage  system,  as  shown 
in  Figure  1.  The  elapsed  time  from  receipt  by  the  Branch  to  the 
start  of  processing,  and  from  receipt  by  the  Branch  to  completion  by 
the  analyst,  constitutes  the  waiting  time  and  turnaround  time,  respec- 
tively. The  service  time  consists  of  the  actual  time  the  analyst 
spends  processing  an  individual  case. 


3.  Applicable  Model 

In  the  model  described  in  this  paper  for  a multiple - 
station,  single-stage  system,  each  station,  k,  is  an  individual 
analyst  where  k>l,  and  each  case  is  processed  by  a single  station 
or  analyst.  Queueing  models  for  this  type  of  system  involve  consi- 
deration of  the  mean  arrival  rate  (\)  and  the  mean  service  or 
processing  rate  to  determine  the  expected  waiting  time  which, 

when  added  to  the  mean  sei*vice  time,  provides  the  total  expected 
time  in  the  system.  The  model  assumes: 

a.  Arrivals  are  random  and  are  Poisson  distributed  with  a 
mean  number  of  arrivals,  \ , per  unit  of  time. 


b.  Service  times  are  constant  or  exponentially  distributed 
with  a mean  service  time,  1//J,  ; i.e.,  the  reciprocal  of  the  mean 
service  rate. 


c. 


The  sei*vice  priority  is  first  come,  first  sei*ved. 


d.  The  number  of  service  stations  (analysts)  is  known,  and 
all  stations  have  the  same  service  capacity. 


e.  The  mean  service  time  (processing  speed)  is  unaffected 
by  the  length  of  the  queue;  i.e.,  service  times  are  statistically 
independent. 


f.  There  is  an  infinite  source  of  inputs. 


Given  that  these  assunptions  are  true,  the  queueing  models  will 
provide,  on  a probabilistic  basis,  measures  of  the  system  character- 
istics that  are  suitable  for  management  analysis  and  decision. 


The  utilization  parameter  of  the  system  is  defined  as  the  pro- 
portion of  time  that  all  k analysts  are  being  utilized,  and  assumes 
that  k jU,  > X . Note  that  if  this  assuirption  is  incorrect  that  the 
k analysts  will  be  theoretically  busy  more  than  100  percent  of  the 
time,  and  an  infinite  queue  will  build  ip.  The  parameter,  denoted 
by  the  Greek  letter  rho,  is  expressed  as 


= X 

w 
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Q INPUT  (FMS  CASE) 


SERVICE  STATION  (ANALYST) 


FIGURE  1. 


FMS  CASE  PROCESSING  SYSTEM: 


Multiple-Station,  Single-Staqe 
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The  probability  that  all  k analysts  are  idle,  i.e.,  the  probability 
that  there  is  no  case  in  the  system,  is  defined  as 


P 

o 


k-1 
L i=o 


1 

f X 1 i 

4. 

1 

xl 

i ! 

M / 

' 

k ! 

The  probability  that  a randomly  arriving  case  has  to  wait  for 
processing,  i.e.,  the  probability  that  there  are  k or  more  cases  in  the 
system,  is 


P 

n> 


k 


(k-l)  ! (k|U  - A ) 


The  mean  (expected)  length  of  the  queue,  i.e.,  the  number  of  cases 
waiting  for  processing  is 


E(m)  = P^ 


X UCA/U)^ 

(k-l)  ! CkjU-X)^ 


The  mean  (expected)  number  of  cases  in  the  system  is 


E(n)  = E(m)  + 


X 


The  average  (expected)  waiting  time  of  an  arriving  case  in  the 
queue,  expressed  in  the  same  time  units  as  \ and  jU  , is 


E(w)  = P 


U CX/U)^ 


(k-l)  ! (kjU-X)' 


= E(ml 

X 


The  average  (expected)  time  a case  spends  in  the  system,  i.e.,  the 
expected  turnaround  time,  expressed  in  the  same  time  units  as 
X and  jJi  , becomes 

E(t)  = E(w)  + 

4.  Data  Collection  and  Analysis 


To  provide  inputs  to  the  models,  the  arrival  and  processing 
distributions  were  analyzed  as  follows: 


a.  Data  for  numbers  of  arrivals  p>er  day  were  collected  for  the 
period  2 Jan  through  14  Apr  76,  and  are  shown  in  Table  1 by  month. 

The  distribution  of  observed  frequencies  of  arrival  determined  from 
the  data  is  also  shown  in  Table  2.  The  cumulative  mean  arrival  rate 


TABLE  1.  FMS  CASE  ARRIVALS  (2  JAN  - 14  APR  76) 


No.  Cases 

Number  of  Days 

Per  Day 
— g — 

January 

February 

March 

1 

April 

TOTAL 

1 

1 

2 

X 

1 

1 

±. 

4 

2 

2 

3 

2 

7 

3 

5 

3 

1 

1 

10 

4 

8 

5 

1 

1 

15 

5 

4 

1 

4 

1 

10 

6 

3 

5 

2 

10 

7 

2 

4 

6 

8 

2 

4 

1 

7 

9 

1 

1 

10 

1 

1 

11 

12 

13 

14 

1 - 

1 

TOTAL 

21 

19 

23 

10 

73 

SOURCE:  DRSAR-CPE-D  daily  logs,  2 Jan  - 14  Apr  76. 


G 


TABLE  2.  FREQUENCY  DISTRIBUTION  - FMS  CASES  (2  JAN  - 14  APR  76) 


No.  of  Cases 
Arriving  Per  Day 
(X) 


Observed 

Frequency 

(fo) 


Expected 
Frequency  (£  ) 
C X=4.80)® 


0 

1 

0 

0.60 

1 

4 

4 

2.88 

2 

7 

14 

6.92 

3 

10 

30 

11.07 

4 

15 

60 

13.30 

5 

10 

50 

12.75 

6 

10 

60 

10.21 

7 

6 

42 

7.00 

8 

7 

56 

4.20 

9 

1 

9 

2.24 

10 

1 

10 

1.07 

11 

0 

0 

0.47 

12 

0 

0 

0.19 

13 

0 

0 

0.07 

14 

1 

14 

0.02 

15  or  more 

0 

0 

0.01 

TOTAL 

73 

349 

73.00 

Mean  Arrival  Rate 


349 

73 


= 4.78  cases/day 


f = T£ 
e ^ o 


P(r=x|  X “ 4.80) 
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was  then  calculated  as  shown.  Since  the  mean  arrival  rate  parameter, 

and  the  total  frequency,  , are  sufficient  to  completely 

describe  a Poisson  distribution,  expected  frequencies,  f , for  a 
Poisson  distribution  with  X"  4.80  were  calculated  using  a table  of 
Poisson  probabilities  of  the  type  found  in  most  basic  statistics  texts, 
and  are  also  listed  in  Table  2.  The  value  of  4.80  was  the  closest 
table  value  to  the  calculated  value  of  4.78.  Observed  and  expected 
frequencies  were  plotted.  Figure  2,  showing  that  arrivals  over  the 
period  plotted  are  closely  approximated  by  a Poisson  distribution  with 
X=  4.80.  A statistical  test  of  significance,  called  the  chi-square 
test,  can  be  used  to  determine  the  goodness  of  fit  of  the  distribution 
of  observed  frequencies  to  the  Poisson  distribution,  and  is  described 
in  Appendix  A for  the  analysis  of  FMS  case  arrivals. 

b.  Reliable  observed  or  estimated  processing  time  per  case 
data  were  not  available.  It  is  known  that  processing  times  are  not 
uniformly  distributed;  therefore,  an  exponential  distribution  was 
assumed.  Since  arrivals  fit  the  model  assumption,  this  assumption 
of  service  time  distribution  was  considered  reasonable.  The  single 
parameter  required  to  describe  an  exponential  distribution  is  the 
mean  processing  time,  1/fJi.  An  estimate  of  jU  was  obtained  by  dividing 
the  total  hours  expended  on  FMS  cases  during  the  period  January  1976 
through  March  1976  by  the  total  number  of  cases  reviewed  during  the 
same  period.  These  data  are  shown  in  Table  3.  Estimates  of  fJ.  are 
also  shown  by  month,  and  cumulatively  through  the  period  analyzed 

(see  paragraph  C4c).  Conversion  to  a mean  processing  rate,  in  cases 
per  day  per  person,  was  necessary  to  provide  input  to  the  model  to 
match  the  mean  arrival  rate  parameter  in  cases  per  day.  The  processing 
times  and  rates  obtained  are  gross  values  in  that  all  times  expended 
in  both  direct  and  indirect  support  of  case  reviews,  including  resolu- 
tion of  policy,  meetings,  and  administrative  activities,  is  included. 

c.  The  arrival  and  processing  data  were  analyzed  to  determine 
the  stability  of  the  derived  parameters  over  the  period  covered  by  the 
data.  Both  monthly  and  cumulative  mean  arrival  and  processing  rates 
were  calculated  as  shown  in  Table  4 and  3,  respectively.  These  analyses 
show  that,  over  the  periods  covered,  the  mean  arrival  rate  increased 
while  the  mean  processing  time  decreased.  These  trends  are  plotted  in 
Figure  3.  It  would  be  expected  that  further  investigation  would  indi- 
cate the  reasons  for  these  trends.  For  exanple,  training  of  analysts 
over  the  period  studied,  with  attendant  learning,  might  explain  the 
downward  trend  of  mean  processing  time. 

5.  Solution  of  the  Model 

The  model  for  multiple- station,  single-stage  waiting  line 
systems  has  been  programmed  on  the  WANG  700  Series  progranriable  calcula- 
tor located  in  the  Cost  Analysis  Division.  The  inputs  are  k,  the  number 
of  stations  (analysts),  in  whole  integers;  mean  arrival 

rate,  in  cases  per  day;  and  jU  (mu) , the  mean  processing  rate  per  person, 
in  cases  per  day.  When  ’kjj.  <.  the  program  will  not  execute. 
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TABLE  3.  FMS  CASE  PROCESSING  HOURS  (JAN  - MAR  76) 


, , January 

Man-Hours  Expended  -J,  586 

No.  Cases  Reviewed  — 68 


February  March  Total 

667  626  1,879 

88  114  270 


Monthly 
Jan  76: 
Feb: 
Mar: 


Mean  Processing  Times  (1/jU ) and  Mean  Rates  (jU) 


586 

68 


= 8.618  hours  per  case  (0.928  cases  per  day  per  person) 


667 

88 


= 7.580  hours  per  case  (1.055  cases  per  day  per  person) 


626 

114 


= 5.491  hours  per  case  (1.457  cases  per  day  per  person) 


Cumulative  Mean  Processing  Times  (l//i)  and  Mean  Rates  (jU) 

Jan  76:  See  Above 
1253 

Jan-Feb;  ~ 8.032  hours  per  case  (0.996  cases  per  day  per  person) 

1879 

Jan-Mar:  ~27o  6.959  hours  per  case  (1.150  cases  per  day  per  person) 

\J  Source:  Labor  Tally  (Detail),  ARICOM  Form  1. 

TJ  Source:  DRSAR-CPE-D  periodic  reports. 
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TABLE  4.  FMS  CASE  MONmY  AND  CUMULATIVE 


MEAN  ARRIVAL  RATES  (2  JAN  - 14  APR  76) 


Month (No.  Days) 

Jan  76  (21) 

Feb  (19) 

Mar  (23) 

1-14  Apr  (10) 


No.  of  Arrivals 

73 

88 

122 

66 

“3^ 


Mean  Arrivals 
Per  Day(X) 

3.48 

4.63 

5.30 

6.60 


Jan  76 
Jan- Feb 
Jan-Mar 
Jan- 14  ^r 


Cumulative  Mean  Arrival  Rates  (\) 
73 

^ * 3.48  cases  per  day 

=»  4.03  cases  per  day 
283 

-gj  * 4.49  cases  per  day 
349 

—sj  = 4. 78  cases  per  day 


Average  No.  of  Arrivals  Per  Oay/Average  Processing  Hours  Per  Case 


FIGURE  3.  FMS  CASE  MEAN  ARRIVAL  RATE  ANO  PROCESSING  TIME 
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HOURS 
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1 
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The  outputs  are  the  systan  characteristics  defined  in  paragraph  C3. 

The  program  facilitates  iterative  solution  of  the  model  for  various 
values  of  the  inputs,  and  is  expecially  useful  viien  trends  exist 
as  discussed  in  paragraph  C4c.  In  this  situation,  management  must 
decide  if  the  system  parameters  have  stabilized  or  if  further  analysis 
is  necessary.  As  an  interim  solution,  the  model  has  been  solved  using 
values  of  \ and  jll  for  both  the  entire  period  covered  by  the  data  and 
for  the  month  of  March  only,  and  increasing  values  of  k.  Printouts 
of  both  solutions  are  included  as  Appendix  B. 

6.  Analysis  of  the  Solution 

Since  the  mean  arrival  rate  is  independent  of  vdiat  happens 
in  the  processing  system  itself,  queueing  analysis  can  only  indicate 
how  waiting  time,  E(w),  and  total  time  in  the  system,  E(t)  , is  affected 
by  changes  in  the  number  of  analysts,  k,  and/or  in  the  mean  processing 
rate,/J..  Given  a stable  value  of  the  latter  as  well  as  the  arrival 
rate,  adding  resources  can  only  reduce  the  total  expected  (turnaround) 
time  to  that  of  the  mean  processing  time,  1//-L.  The  solutions  of 
Appendix  B show  this  to  occur  at  k * 10  and  k = 11,  when  the  input 
parameters  are  for  March  and  the  total  period  studied,  respectively, 
reading  E(t)  to  three  decimal  places.  These  values  are  obviously 
unrealistic,  if  only  because  utilizations  are  less  than  40  percent  at 
these  staffing  levels.  It  is  then  helpful  to  managanent  decision-making 
to  plot  expected  turnaround  time  as  a function  of  k.  These  are  shown 
in  Figure  4.  Both  plots  show  graphically  that  turnaround  time  drops 
rapidly  frcm  the  lowest  value  of  k at  which  kjLX>\as  analysts  are 
added.  Once  management  has  determined  stable  values  of  the  input 
parameters,  this  type  of  analysis  is  very  helpful  in  determining  the 
staffing  level  required  to  achieve  the  desired  turnaround  time,  and 
the  cost  penalties  of  achieving  small  additional  incranents  of  reduced 
time. 


Solution  of  the  model  also  shows  significant  changes  in  turnaround 
time  achievable  by  reduction  in  the  mean  processing  time,  vdien  such 
reductions  are  possible.  For  example,  using  the  March  arrival  rate 
for  FMS  cases,  if  the  mean  processing  rate  can  be  increased  from  1.46 
cases  per  day  per  person  to  2.00  cases  per  day  per  person,  the  calculated 
turnaround  time  at  k = 4 is  reduced  by  71  percent. 


E(t)=Expected  Turnaround  Time,  In  Days 


FIGURE  4 

2.25 
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1.5 


1.0 


0.5 


EXPECTED  TURNAROUND  TIME,  IN  DAY5,  AS  A FUNCTION  OF  THE  MIJMPPR 
OF  ANALYSTS 


NOTE;  1.  Upper  curve  based  on  mean  arrival 

rate,  \ ■ 4.78,  and  mean  processina 
rate,U  = 1.15;  lower  curve  based 
on  A = 5.30  and  jU  = 1.46. 

2.  Lower  limit  is  where  k=5  and  k=4 
for  upper  and  lower  curve,  respec- 
tively, since  queue  length  becomes 
infinite  when  \ > kjJi . 
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APPENDIX  A 


CHI-SQUARE  SIGNIFICANCE  TEST 


CHI-SQUARE  TEST  OF  GOODNESS  OF  FIT 


In  testing  for  goodness  of  fit,  it  is  necessary  to  conpare 
the  distribution  of  observed  samples  with  the  theoretical  distribution 
that  is  assumed  to  be  the  population  from  vdiich  the  sample  was  drawn. 

If  there  is  a high  degree  of  conformity,  any  small  differences  between 
the  observed  and  theoretical  distributions  are  assumed  to  be  the  result 
of  sanpling  variation.  If  there  are  large  differences,  the  conclusion 
may  be  made  that  the  observed  sanple  data  was  drawn  from  a population 
distribution  other  than  that  assumed  in  the  test. 

In  the  case  of  queueing  models,  arrivals  are  assumed  to  be  Poisson 
distributed.  For  FMS  cases,  goodness  of  fit  of  the  observed  arrival 
data  to  the  theoretical  distribution  was  determined  graphically  by  the 
plot  of  Figure  2.  However,  if  it  is  desired  to  test  the  goodness  of 
fit  statistically,  the  chi-square  test  may  be  used.  The  chi-square 
distribution  is  expressed  as 


X 


2 


2 


where  f = an  observed  frequency,  and  f = a theoretical  (expected) 
frequenSy.  The  distribution  characteriltics  are  completely  defined 
by  the  number  of  degrees  of  freedom  (d. f . ) , which  can  be  defined  as 
the  number  of  frequency  classes,  or  groipings  of  data,  being  compared, 
minus  one  d.f.  for  each  restriction  placed  on  the  expected  distribution. 
In  the  case  of  the  Poisson  distribution,  since  it  can  be  conpletely 
determined  by  \ and  the  total  frequency,  there  will  be  two  restrictions. 


Four  conditions  should  be  met  in  order  to  apply  a valid  chi-square 
test.  The  first  two  state  that  the  sample  observations  should  be  inde- 
pendent of  one  another  and  drawn  frcm  the  population  being  analyzed, 
and  that  the  data  are  usually  of  some  nominal,  or  moderate,  measurement. 
The  remaining  two  place  seme  quantitative  restrictions  on  the  data  and 
resulting  frequency  classes,  in  that  the  sample  should  contain  no  fewer 
than  fifty  observations  and  groiped  with  at  least  five  observations  in 
each  frequency  class.  Statisticians  differ  on  the  exact  numbers  in  the 
latter  two  conditions.  However,  fifty  or  more  observations  should  be 
possible  for  most  waiting  line  systans,  and  groiping  of  observed  data 
into  classes  of  five  or  more  should  normally  present  no  problem.  Infor- 
mation obtained  by  discussion  with  Mr.  Harold  Gehle  of  the  US  Army 
Management  Engineering  Training  Agency  indicates  that  the  test  may  be 
considered  valid  with  as  few  as  three  observations  in  any  expected  class. 


The  chi-square  test,  as  with  other  tests  of  significance,  is  set 
up  by  establishing  a null  hypothesis,  H , that  the  sanple  is  drawn 
from  the  theoretical  population  distribStion,  and  an  alternate  hypo- 
thesis, H , that  it  is  not.  The  value  of  chi-square  is  calculated 


17 


from  the  sample  data  using  the  formula  on  the  previous  page.  A test 
criterion  is  then  established  upon  which  to  base  acceptance  of  either 
hypothesis.  The  criterion  is  based  on  the  value  taken  from  a table 
of  values  of  chi-square  of  the  type  found  in  most  basic  statistics 
texts.  The  table  value  is  selected  on  the  basis  of  both  the  number 
of  degrees  of  freedom  and  the  degree  of  risk,  called  Q , that  manage- 
ment is  willing  to  take  that  the  test  statistic  will  exceed  the  table 
value  strictly  by  chance,  thus  leading  to  a rejection  of  (acceptance 

of  H^)  when  it  is  true.  If  is  not  rejected,  we  can  conclude  that 

the  difference  between  the  observed  and  theoretical  distributions  is 
due  to  sampling  error,  since  we  did  not  observe  the  population. 

The  test  procedure  is  best  illustrated  by  application  to  the 
FMS  case  distribution,  at  an  OL  of  0.05: 

1.  H^:  The  population  distribution  of  arriving  cases  is 

Poisson  distributed. 

2.  H : The  population  distribution  of  arriving  cases  is  not 

Poisson  distributed. 

3.  The  number  of  degrees  of  freedom  is  the  number  of  frequency 

classes  minus  the  number  of  requirements  of  the  Poisson 
distribution  ( \ and  total  frequency  = 2) . In  order  to 
meet  the  condition  of  five  or  more  observations  per  class, 
the  data  of  Table  2 must  be  reorganized  as  shown  in  Table 
A. 1.  The  modified  number  of  frequency  classes  is  eight, 
resulting  in  d.f.  = 8-2  = 6.  The  table  value  for  q 05 

(6  d.f.)  = 12.592. 

4.  Criterion:  Reject  H (accept  H ) if  X >12.592;  do  not 

2 ^ ^2 

reject  if  X — 12.592,  when  X computed  using  the 
previously- defined  formula. 

5.  Since  X^  (2.087) < X^  (12.592),  is  not  rejected;  the 
observed  frequency  distribution  is  Poisson. 
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TABLE  A-1. 


FREQUENCY  DISTRIBUTION-FMS  CASE  ARRIVALS 


No.  o£  Cases 
Arriving  Per  Day 
CX) 

£ 

0 

£ 

e 

1 or  fewer 

5 

3.48 

2 

7 

6.92 

3 

10 

11.07 

4 

15 

13.30 

5 

10 

12.75 

6 

10 

10.21 

7 

6 

7.00 

8 or  more 

10 

8.27 

73 

73.00 

(£  -£  ) 
^ 0 e“^ 

(£  -£ 

^ 0 e'^ 

— T" 

e 

1.52 

2.3104 

0.664 

0.08 

0.0064 

0.001 

-1.07 

1.1449 

0.103 

1.70 

2.8900 

0.217 

-2.75 

7.5625 

0.593 

-0.21 

0.0441 

0.004 

-1.00 

1.0000 

0.143 

1.73 

2.9929  2 

0.362 
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APPENDIX  B 


QUEUEING  MODEL  SOLUTIONS 
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□ wanG  LABORATORIES.  INC  6)5  0204  □ WANG  I ABORATC^^INC  6)5  020 


A.l.  SOLUTION  BASED  ON  MARCH  76  DATA 

VbhllLlUG  MObhL  - hbL'i’IPLL  STATION 


K 

- 

4.000 

Xaubua 

Bt 

5.  300 

KUU 

» 

1.460 

rho 

" 

.907 

pO 

. 010262 

Pn,  k 

- 

.603047 

E Lw) 

m 

1.467124 

E (m) 

- 

7 .661761 

E (u) 

* 

11.511696 

ii(t) 

as 

2.172056 

liULULIWG  HUUEL  - ilULTIPLL  STATION 


a 

iauuua 

mu 

rhu 


000 
5.  300 
i.  AfaO 
. 726 


pO 

Pn,K 
iiLw) 
£,  (,ia) 
ji  (n) 
h(t) 


.021926 
.420474 
.210237 
1.114257 
4. 744394 
.695166 


I 


OULOElwO  MODjiL  - llbLlIPLL  STATION 


K.  “ 

6.000 

lambda  ■ 

5.300 

W 

mu  ■ 

1.  460 

r bo  “ 

. 605 

P“  « 

. 025171 

P n , k “ 

.202552 

ii  (v ) “ 

.056541 

iL(m)  - 

.310267 

w 

E (n)  “ 

3.940404 

ii  (,  t)  “ 

. 743472 

OULbElNG  hODLL  - llbLlIPLE  STATION 


k 

M 

7.000 

laubda 

- 

5.300 

lUU 

- 

1.460 

rao 

“ 

.516 

^0 

m 

.026121 

i?a>  k 

- 

. 06  9^*34 

ii(v) 

. 016177 

ii  (m) 

.090342 

(,u) 

n 

3 .72  647  9 

liU) 

- 

. 7031o9 

t 
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i^bLlLlliO  liULiLL  - libLlIPLii  bTATIOb 


i ^ ~ 

) laubda  ■>. 

( mu 

j , rho  -! 

i 

I Pu,k  “ 

bi(w)  - 

. EXm) 

»( 


t-  '• 


laubda  “I 
mu'  “I 
riio^  - 

- - -e  -i ♦ 


- ( 


jjO- 

Pn,k. 
i.XwJ.r- 
h\u)  -•; 
bCn) 
btt) 


d . uuo 
5.3  00 
- 1.4  60 
.433 

. .026^01 
! .030130 
.005666 
! .030030 
-33  r6601-67- 
* . 69  0397 


liUJiUEING  MODEb  - llbLTIPLL 


^ ■!  ‘ 
^ ■ T 


? j 


y^aoo 

5.300 
1.400 
.403 


1 .... 


f • 


• I 

i • 

f 

T i 

3-i: 
1 • 
I 


I 


.026462 
.013369 
.001707. 
T .009051 
3 . 639166 
I .666639 
4. 


K 

laubda 

lilU 

r hu 


I 

i 

liUEliElIvG  hOOLL,  - labLIIPlL 

' : i ' ' ’ ; 

10.000 
5..  300 
1.460- 
i .363 


Pn,  k 
6 (.w) 
L (m) 
li  (u) 
Lit) 


.026304 

.004356 

.000469 

.002396 

.632733 

.665421 


v^biibLlbG  liOULL  - libLlll’LL 


K. 

= ‘ 

11.000 

lauDua 

a 

3.300 

uu 

1.460 

r tio 

sa 

. 330 

po 

ss 

.026j10 

Pq,  K 

St 

.001430 

li  (.v) 

.000132 

Mm) 

» 

.000704 

Liu) 

S3 

3 .030641 

hit) 

XX 

. 66  5 0 64 

6 TAT 10 b 


STATION 


STAllUlJ 
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WANG  LABOR ATORIES.  INC  615  0204  O wanG  LABORATORIES  INC  6I5-0204 


tilLbLlIlG  HOULL  - HULTITIL  tiTAlIUN 


J 

K ” 

12.000 

laubila  = 

5. 300 

JilU  •• 

1.460 

r ho  ■ 

.302 

. 02d312 

Pn , K ■ 

. 000415 

h (v)  “ 

. 000034 

L.(iu)  «= 

. 0 OCXu  0 

mil)  = 

3.630317 

mt)  “ 

. 664965 

OOLOLILG  riu] 

K “ 

13 .000 

w 

iaiabaa  *» 

J.  300 

mu  “ 

1.460 

riiO  « 

.279 

m 

. 026  j 12 

•= 

. 000112 

li  (.w)  = 

. OOOOOo 

M.ia)  - 

. 000043 

hUi)  “ 

3.63  Olo  0 

iiLt)  - 

. 6o4939 

IILLIII’LL 


l^llLbLlbb  IlObi^L  - ilbLlirLL 


S 

14.000 

laubua 

KS 

5.300 

mu 

B 

1.460 

r ao 

BS 

. 25  y 

po 

B 

.026512 

P n , k 

S 

.000020 

'W' 

h (,v) 

B 

.000001 

L (u) 

BS 

. 000009 

£.  Cll) 

« 

3 .630146 

E(t) 

m 

.684933 

I 

i 


l^ULbAlbG  hUbLL  - MbLlirLL 


k - 

iauDcia  ■ 
luu  •• 
rno  ” 


li. 000 
5.300 
1.460 
.242 


a 


pO 

I'n,  K 
ii(w) 
n (w; 
h (u) 

i^Lt) 


.02OJ12 
. 000006 
. 000000 
. 000002 
3 . 63013  9 
. 604931 


28 


station; 


SlATIUlt 


STATION 


t 


l^ULLLIUG  llUDLl. 


MLLTIPLL  STATION 


k *; 
luutda  «i 
mu  " 
rho 

I 

i p o 

Pn,  k ■ 

• - fcr(w>-  -j 
, li'lto)  “i 
iii(n) 
k'(,t)  -I 


4.  UUO 

5.  300 
2.000 
, .b6  2 

• .061137 
. ..372340 
f .1*37  906- 
! . 730903 
3.300903 
t .637906 


I 

I 
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A. 2.  SOLUTION  BASED  ON  2 JAN- 14  APR  76  DATA 


i^ULbLll.G  HOLLL 


iiLLiirLL  bii/noi; 


iT> 

lU  ■'W' 

k 

« 

0.  uuu 

QC 

o 

laidbua 

e 

4. 7bC 

< 

tuu 

sai 

1.  ibO 

o 

< ^ 

r no 

“ 

.b3i 

c 

2 

P* 

.UlOUbl 

< 

Pu,  k 

B 

. 01 6u/7 

o 

b (.w) 

= 

.03307  9 

L (to) 

B 

3. 035O02 

(u) 

7.192204 

k (t) 

■ 

1. 3 U4  043 

qUiiULliMG  hUULL  - nLLllPLii  bXAlIUb 


k - 

0. 000 

iaiiibaa  “ 

4. 7o0 

W 

mu  » 

1.150 

r ho  “ 

.092 

f B 

.013917 

P u , k “ 

. 324440 

0(W)  “ 

.133037 

k (to)  = 

.731320 

W 

k(n)  - 

4.b0bo42 

k(t)  “ 

1.022003 

i^biibLllMG  MUbLL 


libLTIPLL  blATIUb 


O 

VO 


k 

« 

7 . 000 

iautda 

- 

4.700 

uu 

1.150 

r iio 

" 

.393 

po 

- 

.013112 

Pii,k 

. 1 3 u 2 2 0 

n (v) 

m 

.040303 

k (to) 

- 

.231202 

L (u) 

« 

4.307004 

k (t) 

“ 

.9179jO 

C^UaLLIInG  1-lObLL  - ULLTIl’LL  STAlIUb 


K. 

0.000 

laubaa 

“ 

4.  700 

mu 

1.150 

r hu 

.319 

pp 

B 

.013491 

i:^n,  k 

B 

. 071247 

k(w) 

as 

.010119 

k (to) 

SB 

.07  7030 

k (u) 

» 

4.233571 

k(t) 

“ 

.005004 
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1<ULLE1NG  hUDLL  - HbLTlTLL  STaTIOII 
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sx 
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lambda 

■1 

4. 7bU 

mu 

M 

1.150 

rho 
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P* 
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ii  (.a) 

* 
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b(t.) 

: .074079 

qUliVKlMG 


hUULL 

i. 


- IILLIIPLE  STATIOI, 


k. 

mr 

j 

10.000 

lambda 
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4.700 

mu 

* 

il.150 
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P* 
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B 

’ .002450 
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- 

4.150971 

k(t) 

m 

, .070077 

- lILLTirLL  STATION 
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k 
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QLLLLiUO  MODEL  - MLLTIPLL  D1ATI01, 
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Tlication  of  queueing,  or  waiting-line,  theory  to  the  the  review  and  validation 
^ign  Military  Sales  (FMS)  cases  by  the  Cost  Analysis  office  of  an  Army  Materiel 
iment  and  Readiness  Command  major  subordinate  command  is  presented  in  this  report 
'>q  analysis  is  applied  to  determine  the  behavior  characteristics  of  the  pro- 
) system  so  that  management  can  take  appropriate  action  to  reduce  the  total  time 
system.  The  primary  characteristics  analyzed  are  the  expected  length  of  the 
the  expected  number  of  cases  in  the  system,  the  expected  waiting  time,  and  the 
-?d  turnaround  time.  Queueing  theory  is  described,  and  the  application  to'  FMS 
^resented,  to  also  serve  as  an  aid  in  the  solution  of  similar  types  of  waitinq- 
^-oblems. 
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